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Abstract 

We use optical tweezers to perform stretching experiments on DNA molecules when interacting 
with the drugs daunomycin and ethidium bromide, which intercalate the DNA molecule. These 
experiments are performed in the low-force regime from zero up to 2 pN. Our results show that 
the persistence length of the DNA-drug complexes increases strongly as the drug concentration 
increases up to some critical value. Above this critical value, the persistence length decreases 
abruptly and remains practically constant for larger drug concentrations. The contour length of 
the molecules increases monotonically and saturates as drugs concentration increases. Measured in- 
tercalants critical concentrations for the persistence length transition coincide with reported values 
for the helix-coil transition of DNA-drug complexes, obtained from sedimentation experiments. 
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I. INTRODUCTION 



DNA-drug interactions have been much studied along the past years. An important 
motivation for these studies is the fact that many of the studied drugs are used for treatment 
of human diseases, particularly, in cancer chemotherapy. 

Single molecule stretching experiments using optical tweezers have yielded agreat amount 
of information about DNA-protein and DNA-drug interactions 1, 2, 3), 4, [J 

Recently, we studied the interaction between psoralen and DNA when illuminated with 
ultraviolet light A (UVA). Psoralen is a drug used to treat some skin diseases, like psoriasis 
and vitiligo. This drug intercalates the DNA molecule and can form covalent linkages with 
the thymines if the complex is illuminated with ultraviolet light, modifying drastically its 
elasticity and impeding the DNA replication and transcription. The persistence length of 
the DNA-psoralen complexes formed after UVA illumination were measured in [3]. 

Daunomycin and ethidium bromide (EtBr) are other examples of drugs which intercalate 
the DNA molecule and can modify its elasticity, dependingon the drug concentration. Both 
drugs unwind the DNA double helix when intercalating [8j. Daunomycin is an anthracy- 
cline antibiotic used in the treatment of various cancers. It inhibits DNA replication and 
transcription when intercalating, impeding cell duplication 9j. Ethidium bromide (EtBr) 
is commonly used as a non-radioactive marker for identifying and visualizing nucleic acid 
bands in electrophoresis and in other methods of nucleic acid separation. 

Several works have reported different results for the effects of these drugs on the entropic 
elasticity of DNA molecules. In those works, the measured parameter used to study elasticity 
modifications is the persistence length of the DNA-drug complex. Smith et al. {2] report 
that ethidium bromide does not modify the elasticity of the DNA molecule, only increasing 
its contour length by ~40%. Tessmer et al. [h]] report that ethidium bromide causes a large 
increase in the contour length and a decrease in the persistence length of the complex for 
1 /j,M of the drug, and at lower concentrations, an increase in both persistence and contour 
lengths. Recently, Sischka et al. H[ report the value of 28.1 nm for the persistence length of 
DNA-daunomycin complexes and 20.7 nm for DNA-EtBr complexes, smaller than the bare 
DNA persistence length of about 50 nm. The authors have used in this work a concentration 
of 1 fiM for both drugs, and a DNA concentration of 15 pM. In the present work, in order to 
clearly establish the effect of these intercalating drugs on the persistence length of the DNA 
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complexes, we performed stretching experiments at various drugs concentrations, from zero 
up to saturation of the complexes. We show that the values obtained for the persistence 
length depend strongly on the concentration ratio between drug and DNA base pairs. Our 
results show that the persistence length of the complexes increases as we increase the drug 
concentration until certain critical concentration is reached. Above this critical concentration 
the persistence length decreases abruptly and remains practically constant for larger drug 
concentrations. 



II. EXPERIMENTAL PROCEDURE 



To measure the persistence and contour length of DNA molecules and DNA-drug corn- 
ilexes, we use optical tweezers [l[ [2I, and intensity autocorrelation spectroscopy 



121. 



The samples consist of A-DNA molecules in a PBS pH 7.4 with [NaCl] = 140 mM solution. 
We attach one end of the molecule to a microscope coverslip, and the other end is attached 

n 

to a polystyrene bead. To do this, we use the procedure described in [13J . We add the drug 
in the sample immediately before the measurements. The DNA concentration used in all 
experiments was Cdna — 6.81 /ig/mL, which corresponds to a base pairs concentration of 
Cbp = 11 fiM. 

Our optical tweezers is mounted in a Nikon TE300 microscope with an infinite corrected 
objective (100X, N.A. = 1.4). The trapping laser is an infrared (IR) laser with A = 832 nm 
(SDL, 5422-H1). The optical tweezers is used to trap the polystyrene bead attached to the 
end of the DNA molecule, so we can manipulate and stretch the DNA molecule. 

In addition, we use a He-Ne laser (A = 632.8 nm) as the scattering probe. The backscat- 
tered light by the polystyrene bead is collected by a photodetector, which delivers pulses 
to a digital correlator. We, then, obtain the autocorrelation function of the backscattered 
light, which allows us to determine the stiffness of the optical trap, due to the Brownian 
motion of the trapped bead. 

The next step in the experimental procedure is to obtain the force versus extension curves 
for the DNA molecules and DNA-drug complexes. To do this, we use the optical tweezers 
to trap the bead with the DNA while pulling the microscope slide, stretching the DNA. 
The backscattered light is collected while stretching the DNA. From the backscattered light 



3 



intensity one obtains the displacement of the trapped bead in relation to its equilibrium 
position, and by multiplying it times the tweezers' stiffness, the force exerted by the DNA 
molecule while it is stretched is obtained. The details about our experimental setup and 
experimental procedure can be found in j?! \l2\. 

Finally, with the force versus extension curves, we use the approximate expression derived 
by Marko and Siggia (Eq. [1]) [ijj] to obtain the persistence and contour length of the DNA 
molecules and DNA-drug complexes, 
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where z is the DNA molecule end-to-end distance, ks is the Boltzmann constant, T is the 
absolute temperature, A is the DNA persistence length and L is the DNA contour length. 

Figure [1] is a typical force versus extension curve obtained with this procedure, for a 
drug-free DNA molecule. Fitting this curve with Eq. [T], we extract the persistence length 
(A = 50 ± 3 nm) and the contour length (L = 16.5 ± 1 fxm) for the A-DNA. These values 



correspond to well-known values reported in the literature [3 



m 



16J. 



III. RESULTS AND DISCUSSION 



In this section we show the results obtained for the two drugs used: daunomycin and 
ethidium bromide. 



A. Daunomycin 

We have performed experiments with DNA-daunomycin complexes for several drug con- 
centrations. In Fig. [2] we show the persistence length (A) of the complexes as a function 
of total daunomycin concentration (Co) for fixed DNA base pairs concentration of C& p = 
11 //M. We denote by Co the total daunomycin concentration used to prepare the sample, 
which is the sum of both the bounded to DNA and the free drug concentration in solution. 

The point which the drug concentration is zero corresponds to the drug-free DNA situa- 
tion with A = 50 nm. 

The behavior of the persistence length A as a function of daunomycin concentration Co 
can be described as follows: it initially increases with Co until it reaches a maximum value 
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(~ 280 nm) at the critical concentration C c £ %ttcal = 18.3 /iM. Then, it decays abruptly to 
around 75 nm and remains practically constant at this value even if we continue to increase 
Co- The contour length increases monotonically from 16.5 ± 1 fim up to the saturation 
value of 21 ± 1 ^m. These mean values are obtained performing an average over many 
different DNA molecules and DNA-drug complexes. Such distribution of contour length 
values was also observed by Mihailovic et al. [17 ]. 

Also, we can estimate the exclusion parameter n (number of total base pairs divided by 
the number of total intercalated drug molecules) from our experimental data. The average 
value of the contour length for DNA-daunomycin complexes obtained when using a saturated 
concentration of the drug increases about 27% relative to drug-free DNA contour length 
(16.5 /mi). This means that when all possible drug molecules are intercalated, the DNA 
increases its contour length by 4.5 /im. Knowing that each intercalated daunomycin molecule 
increases the contour length of the complex by 0.31 nm [8j, we determine the total number of 
intercalated daunomycin molecules, which is around 14500. Finally, the exclusion parameter 
can be obtained by dividing the number of base pairs of the A-DNA (48500) by the number 
of total intercalated drug molecules (14500). We obtain n = 3.3, in good agreement with 
the value 3.5 reported in 9|. 

For comparison purposes, Fig. [3] shows two force versus extension curves (normalized by 
the contour length) for two daunomycin concentrations, before and after the transition. The 
data points in this figure are smoothed for better visualization, i. e., Brownian fluctuations 
are averaged out. 



B. Ethidium Bromide (EtBr) 

The behavior of the persistence length as a function of the drug concentration for DNA- 
EtBr complexes is very similar to the DNA-daunomycin complexes. The difference is that in 
this case the transition occurs at a lower drug concentration (see Fig. H]) for the same DNA 
base pairs concentration Cb p = 11 /iM. The maximum value measured for the persistence 
length of DNA- EtBr complexes is ~ 150 nm, at the critical concentration C^ ttlcal = 3.1 /xM. 
The contour length increases monotonically from 16.5 ± 1 /mi up to the saturation value 
of 23 ± 1 /mi. Again, these values for the contour lengths are averages over many different 
molecules. 
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Repeating the same calculation for the exclusion parameter of EtBr, which increases 
the DNA contour length by 0.34 nm per intercalated molecule 111 ], we obtain n = 2.5, in 
reasonable agreement with the value 2.01 reported in 181 ]. 



C. Equilibrium binding constants 

In our experiments we control the total drug concentration Ct and the total concentration 
of DNA base pairs Cb p - To discuss the elastic properties of the DNA complex formed the 
important parameter to consider is the ratio r between the concentration of bounded drug 
(Cb) per concentration of DNA base pairs (C&, p ) . In order to obtain r, the binding of molecules 



to DNA is analyzed using the neighbor exclusion model [9|. A closed form for this model 
was obtained by McGhee and von Hippel [19J and can be expressed by the equation 



r 

— = Ki(l - nr) 



-I n-l 

nr 



l-(n-l)r ^ 
where r is ratio between the concentration of bounded drug (Cb) per concentration of DNA 

base pairs (Cbp), C/ is the concentration of free drug (not bounded), Ki is the intrinsic bind- 
ing constant and n is the exclusion parameter in base pairs. For a more detailed discussion 
about the neighbor exclusion model, see [lij. 

The concentration of free drug (Cf) can be simply related with the concentration of 
bounded drug (Cb) and the total drug concentration (Ct) through the equation 



C T = C f + C b . (3) 



Using Eq. El and [3] with 
binding constant reported in 



;he determined exclusion parameter (n = 3.3), the intrinsic 
9j for daunomycin, Ki = 7xl0 5 M _1 , the critical daunomycin 
concentration measured in this work (Cp lUcal = 18.3 /uM), and the concentration of DNA 
base pair used in our experiments (Cb p = 11 /jM), we can determine the critical ratio r c , 
which we define as the value of r at the abrupt transition for the value of the persistence 
length. We then obtain the value r c = 0.248. 



Similarly, for EtBr, we use the parameters n = 2.5, Ki = 1.5 xlO 5 M -1 |20j, and C% %tlcal 
= 3.1 fiM determined again from the abrupt change in persistence length. We obtain r c = 
0.131. In Section [HI Dl we compare the values obtained for r c with those reported in the 
literature for a sedimentation experiment. 



It is important to mention that K{ varies with the ionic strength of the solution. The 
values used here are the values for the ionic concentrations used in our experiments. 



D. Interpretation of the DNA-drug complexes elasticity results 

For low drug concentrations, drug intercalation in the DNA molecule increases the rigi dity 



of the complex (see Figs. [2]andH]). This is consistent with the results of Vladescu et al. 21], 
which shows that EtBr stabilizes the DNA double-helix for low drug concentrations. They 
have performed melting experiments with various EtBr concentrations, from zero to 2.5 fiM, 
showing that EtBr intercalation stabilizes the DNA double-helix in this concentration range. 
Therefore, we expect an increase of the persistence length of DNA-drug complexes in this 
low concentration range. Figure H] shows this increase for EtBr, and Fig. [2] shows a similar 
result for daunomycin. 

For high drug concentrations, i. e., above the critical concentration (peak of Figs. [2] and 
Hj), the persistence length of the complexes decays abruptly and remainspractically constant. 

It is well-known that intercalation unwinds the DNA double-helix [8j. Due to unwinding 
and above some drug critical concentration, the complexes can have a helix-coil transition, 
which can cause DNA denaturing as the DNA is stretched, decreasing the persistence length 
of DNA-drug complexes as seen in Figs. [2]and|H The unwinding angle per intercalated EtBr 
drug molecule is approximately 1.7 times greater than that for daunomycin intercalation {sj]. 
Therefore we expect that the transition occurs for EtBr at a lower drug concentration as 
compared with daunomycin, if the same DNA concentration is used. This is confirmed 
experimentally in our data of Figs. [2] and HI 

Sedimentation experiments performed with circular DNA as a function of daunomycin 
and ethidium bromide concentrations display a minimum in the sedimentation coefficient 
S20 at r c = 0.192, for daunomycin and r c = 0.114 for ethidium bromide jsj]. The minimum in 
the sedimentation coefficient S20 is associated with a helix-coil transition, due to unwinding 
of the DNA double- helix by the intercalating drugs 8J. These numbers agree within 15 to 
30% with the values of r c determined from our DNA persistence length measurements. This 
indicates that the abrupt change of the DNA persistence length for both drugs might be 
also caused by a helix-coil transition due to the unwinding of the DNA double-helix as the 
drugs intercalate into it. The reasonable agreement between the critical ratios r c obtained 
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from sedimentation experiments [8| and from our measurements of the persistence length 
transition provides an evidence that a helix-coil transition is probably what we are observing 
in our experiments. 

In addition, it is known that EtBr (and also most intercalating drugs) exhibits multi- 



modality at their interaction with DNA 22j, |23|]. The type of interaction varies with the 



drug concentration. The abrupt transition shown in Figs. [2] and |4] might as well be caused 
by different ways of drug binding to DNA. 



IV. CONCLUSION 



We have made systematic measurements of the entropic elasticity variation of a A-DNA 
molecule when interacting with two drugs, daunomycin and ethidium bromide, as a function 
of their concentrations. Our results show that the persistence length of the DNA-drug 
complexes increases strongly as the drug concentration increases, for low concentrations. 
Above certain critical drug concentration the persistence length decreases abruptly and 
remains practically constant for high drug concentrations. This behavior is quite similar 
for both daunomycin and EtBr, as shown in Figs. [2] and HI Our results suggests that the 
abrupt transition observed in the persistence length might be due to a helix-coil transition 
and denaturing of DNA-drug complexes above the critical concentration, resulting in a 
decrease of the persistence length. 
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Figure Legends 

Figure U\ 

Force as a function of extension for a drug-free DNA molecule. By fitting this curve with 
Eq. [1], we determine the persistence length A = 50 ± 3 nm and the contour length L = 16.5 
± 1 fxm. 

Figure [H 

Persistence length A of DNA-daunomycin complexes as a function of drug concentration 
for fixed DNA concentration (Cb p = 11 fiM). A initially increases with Co until it reaches 
a maximum value (~ 280 nm) at the critical concentration C c ^ lttcal = 18.3 /iM. Then, the 
persistence length decays abruptly to around 75 nm and remains practically constant at this 
value even if we continue to increase the drug concentration. 

Figure 

Force versus extension curves (normalized by the contour length) for two daunomycin 
concentrations. The data Brownian fluctuations are averaged out for better visualization. 
Circles: Co = 20.1 fiM (above the critical concentration) and A ~ 61 nm; triangles: Co = 
18.3 fiM and A ~ 263 nm. Dashed lines are fittings using Eq. [TJ 

Figure [^} 

Persistence length of DNA-EtBr complexes as a function of drug concentration for fixed 
DNA concentration (Cb p = 11 fjM). Here, the transition occurs at a lower drug concentration. 
The maximum value measured for the persistence length of DNA-EtBr complexes is ~ 150 
nm, at the critical concentration C c ^ ttlcal = 3.1 /iM. 
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FIG. 1: 



12 









































I 

f 






















- .* 




- - i 


\ 




- r 


h 













10 



20 



30 



40 



50 



60 



C d OiM) 



FIG. 2: 



13 




14 



50 



- I 

I 

J 


! 

| | 

i i 

! ! 


! ! 

! ! 

; ; 
i i 

i i 

: : 


i 


» ? - 

i i 

! ! 




i i 

i i 



5 10 15 20 

C E (fiM) 



FIG. 4: 



15 



